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Abstract 



The relativistic two-body potentials of constraint theory for systems composed of two 
spin-0 or two spin-| particles are calculated, in perturbation theory, by means of the 
Lippmann— Schwinger type equation that relates them to the scattering amplitude. The 
cases of scalar and vector interactions with massless photons are considered. The two- 
photon exchange contributions, calculated with covariant propagators, are globally free 
of spurious infra-red singularities and produce at leading order 0(q;^) effects that can be 
represented in three-dimensional x-space by local potentials proportional to (a/r)^. An 
approximation scheme, that adapts the eikonal approximation to the bound state problem, 
is deviced and applied to the evaluation of leading terms of higher order diagrams. Leading 
contributions of n-photon exchange diagrams produce terms proportional to (a/r)". The 
series of leading contributions are summed. The resulting potentials are functions, in the 
cm. frame, of r and of the total energy. Their forms are compatible with Todorov's 
minimal substitution rules proposed in the quasipotential approach. 
PACS numbers: 03.65.Pm, ll.lO.St, 12.20.Ds, ll.80.Fv. 
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1 Introduction 



The three-dimensional reduction of the Bethe— Salpeter equation [|T], ^ was achieved in 
the past with the techniques of the quasipotential approach [||, ^ ^ ^ |, [1^, [11], |12[. The 
method consisted in iterating the Bethe— Salpeter equation around a three-dimensional 
hypersurface in relative momentum space; the result is a three-dimensional wave equation 
for the relative motion, where the potential is related to the scattering amplitude by 
means of a Lippmann-Schwinger type equation. This wave equation is not unique in 
form and depends among others on the three-dimensional hypersurface chosen for the 
projection operation, reflecting the way of eliminating the relative energy and relative 
time variables from the initial four- dimensional equation. A variant of this approach 
consisted in iterating the Bethe— Salpeter kernel around a simpler three-dimensional one 

Quantitative applications of these three-dimensional methods concerned mainly QED 
and the positronium system. It turns out that the Coulomb gauge is the most convenient 
gauge to achieve the perturbative calculations JlT], [l^, |14|, |l6l. Covariant gauges 



T7j, or covariant propagators for scalar exchanges |T^, have the tendancy to produce 
spurious infra-red singularities that are only cancelled by contributions of higher order 
diagrams. This feature considerably reduces the domain of practical applicability of the 
Bethe— Salpeter equation in covariant perturbative calculations. Actually, this flaw can 
be traced back to the noncovariant nature of the three-dimensional equation around which 
iteration is accomplished. Also, the instantaneous approximation [^], applied to covariant 
propagators, does not produce the correct 0{a^) terms in perturbation theory, neither the 
correct infinite mass limit, unless an infinite number of higher order diagrams are taken 
into account [pOf . 

In this respect, the variant of the quasipotential approach, developed by Todorov |I0| , 



where, to some extent, potentials may appear in local form in x-space with appropriate 
cm. total energy dependences, provides a manifestly covariant framework that is free 
of the abovementioned diseases. In particular, it has been shown for one spin-| 



and one spin-0 particle systems, that in the two-photon exchange diagrams, calculated in 
the Feynman gauge, the spurious infra-red singularities cancel out and correct results are 
obtained for the O(a^) and 0{a^\na~^) effects. This approach could not, however, reach 
sufficient generality, because of the lack in it of covariant wave equations for two-spin- 1 



particle systems. 

It is the manifestly covariant formalism of constraint theory applied to two- 
particle systems ^ |25|, that provided the general framework for setting up covariant 
wave equations with the correct number of degrees of freedom. Two-spin--^ fermion sys- 



tems are described there by two compatible Dirac type equations p6|, In local 

approximations of the potentials, these equations, without loss of their covariance prop- 



erty, can be reduced to a single Pauli-Schrodinger type equation |^ which displays their 
practical applicability to spectroscopic problems. The connection of these equations with 
the Bethe— Salpeter equation is obtained with standard iteration techniques [^, |3l|. (For 
a comparison of various three-dimensional equations see Ref. |3^ .) 

The purpose of the present paper is to investigate the perturbation theory properties 
of the potentials of constraint theory wave equations. We shall not consider here radia- 
tive corrections, which actually are calculable independently from any three-dimensional 
formalism, but rather shall analyze the structure of multiphoton-exchange diagrams. 

Our analysis is accomplished in three steps. Firstly, we study the structure of two- 
photon exchange diagrams and show that the corresponding potential is free of spurious 
infra-red logarithmic singularities and yields the correct O(a^) effects; it can be repre- 
sented, in three-dimensional x-space, by means of a local fuction, proportional to (a/r)^. 
Secondly, we device an approximation method for the evaluation of the above diagrams 
that globally reproduces the same result for their leading effect. This method is a variant 



of the eikonal approximation []33| , |34| , p5| , adapted to the bound state problem. 

Thirdly, we apply the approximation method to the higher order diagrams. We show 
that the n-photon exchange diagrams yield for the corresponding potential, as a dominant 
contribution, an 0{a'^"') effect, represented by a term proportional to (a/r)". We then 
sum the series of dominant contributions to obtain the potential in compact form, as a 
function of (a/r) and of the cm. total energy. 

The expression we find for the timelike component of the electromagnetic potential 
can naturally be completed to include the spacelike components as well. It then represents 
the fermionic generalization of the potential proposed by Todorov in his quasipotential 
approach for two-spin-0 particle systems, on the basis of minimal substitution rules |jlO| . 



The same calculations, applied to scalar photons, yield a potential that is equivalent, in 
the static approximation of its higher order terms, to the form proposed by Crater and 
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Van Alstine P^, who extended Todorov's substitution rules to scalar interactions. 

The above calculations result in sizable effects mainly in strong coupling problems. 
One such problem is provided by the evaluation of the short-distance potential in QCD, 
which is very similar to the QED potential, with values of the coupling constant of the 
order of 0.5; furthermore, part of the radiative corrections could be taken into account 
by the replacement of the coupling constant by the effective coupling constant, obtained 
from renormalization group analysis. 

Another domain of application of these potentials is provided by QED in its strong 
coupling regime, where lattice calculations signal a chiral phase transition for values of a of 
the order of 0.3-0.4 Other qualitative results obtained from the Bethe— Salpeter 

equation with one-photon or one-gluon exchanges (such as spontaneous breakdown of 
chiral symmetry) might be tested with these potentials and their stability in the presence 
of multiphoton or multigluon exchanges checked. Also, these potentials could naturally 
be generalized for the treatment of other interactions, such as those corresponding to the 
exchanges of massive particles or to confining interactions, by the replacement there of 
the Coulomb potential by another appropriate potential. 

The paper is organized as follows. In Sec. 2 the wave equations of constraint the- 
ory are introduced. In Sec. 3, their connection with the Bethe— Salpeter equation and 
the relationship of the potential with the scattering amplitude are established. In Sec. 
4, the perturbative calculational method of the potential is developed. In Sec. 5, the 
potential corresponding to the exchanges of two scalar photons for two-spin-0 particle 
systems is calculated to order a^. In Sec. 6, similar calculations are done for scalar and 
vector photons for a fermion-antifermion system. In Sec. 7, the approximation method 
leading to the results of Sees. 5 and 6 is deviced and some of its general consequences 
about cancellation of diagrams are shown. In Sees. 8 and 9, the dominant contribu- 
tions coming from multiphoton exchanges are calculated and summed for two-boson and 
fermion-antifermion systems, respectively. In Sec. 10, the Bethe— Salpeter wave function 
is reconstituted in terms of the constraint theory wave function in the framework of the 
present approximations. Conclusion follows in Sec. 11. Details of the calculations of 
two-photon exchange diagrams are presented in Appendices A and B. A summary of the 



calculation of the vector potential was presented by the present authors in Ref. 
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2 Two-body wave equations of constraint theory 



A system of two-spin-0 particles, with masses mi and m2, respectively, is described by 
means of two wave equations that are generalizations of the Klein— Gordon equation 



23, 24, 25, 27, 28 



{pj-ml - V) ^(xi,X2) = 
{pl~ml-V) ^ixi,X2) = 



(2.1a) 
(2.1b) 



where \^ is a Poincare invariant function or operator of the variables of the system. It is 
possible to choose different potentials in the two equations, but by means of canonical (or 
wave function) transformations one can bring the corresponding equations to the above 
configuration p9| . 

We use the following definitions for the total and relative variables: 



P = Pl+P2 



P = ^{Pi-P2) , M 



mi + 1712 



-^=9(^1 + ^2), X = X1-X2 



(2.2) 



For states that are eigenstates of the total momentum P we define transverse and longi- 
tudinal components of four-vectors with respect to P: 



q.P 
~P^ 



P 



{q.P)P,, P, 



P„ 



q.P 



Pt. 



Ql 



,T2 



q 



-X 



T2 



[2.3) 



This decomposition is manifestly covariant. In the cm. frame the transverse components 
reduce to the three space components, while the longitudinal component reduces to the 
time component of the corresponding four-vector. 

The compatibility condition of the two wave equations ( |2.1a|) -( p.lb|) imposes the fol- 
lowing constraints on the wave function and the potential: 



(pI - pi) -{ml-ml)]^ = , 



pI-pI,v 



. 



The solution of Eq. 



is (for eigenfunctions of the total momentum P): 



-iP.X-i{m\ - ml)xL/{2PL) 



^jj{x'^ 



[2A) 
[2.5) 



(2.6) 



tells us that V does not depend on the relative longitudinal coordinate 

V = \/(x^,Pl,/) . (2.7) 
allow us to eliminate the longitudinal relative coordinate and 



while Eq. 

XL- 



Thus, Eqs. 

momentum. Equation ( |2.4| ) can be rewritten in a slightly different form by defining the 
constraint C{p): 



C{p) = 2PlPl - {ml - m. 



. 



It yields for the individual logitudinal momenta the expressions: 



PlL 



2 



Pl 

V2L = - - 



m\) 



2Pl ' 2 2Pl 

The eigenvalue equation satisfied by the internal wave function i/) is: 



p2 1 

(m, + m,) + 

4 2^^ 2; 



[ml — mlY 



4P2 



+ p 



T2 







:2.8) 



(2.9) 



(2.10) 



which is a three-dimensional Klein— Gordon— Schrodinger type equation. Also notice that 
when constraint C (|2.8|) is utilized, the two individual Klein— Gordon operators become 
equal: 

p2 1 



Ho ^ {pi 



C 



iP2 - ^2, 



C 



(2.11) 



In general the potential V [Eq. ( p.7| )] is a Poincare invariant integral operator in 



X 



but in simplified cases it can be represented as a function of with a general dependence 
on Pl, and an eventual quadratic dependence on p^ . 

The scalar product for the wave functions satisfying Eqs. ( ^.laD -( p.lb| ) can be obtained 
either from the construction of tensor currents satisfying two independent conservation 
laws with respect to Xi and X2, respectively, or directly from the Green's function of 
Eqs. (|2.1a| )- (|2.1b|) or (|2.10|) [§, 0. The result for the norm is (in the cm. frame): 

d 



gp2 



(Ho - V) 



2Pn 



(2.12) 



where Hq is the Klein— Gordon operator (|2.11| ). 

For a system of two-spin- 1 particles, composed of one fermion with mass mi and one 
antifermion with mass m2, the constraint theory wave equations can be written in the 



form 28 



(71 -Pl - "^1)^ = (-72-P2 + m2)V'^ 
[-12-P2 - m2)^ = (71 .pi + mi)V^ , 



(2.13a) 
(2.13b) 



where is a Poincare invariant potential. Here, \E' is a sixteen-component spinor wave 
function of rank two and is represented as a 4 x 4 matrix: 



^aiQ2(a;i,a;2) (ai, as = 1, ■ ■ ■ , 4) 



(2.14) 



where cti (as) is the spinor index of particle 1 (2). 71 is the Dirac matrix 7 acting in the 
subspace of the spinor of particle 1 (index ai); it acts on \E' from the left. 72 is the Dirac 
matrix acting in the subspace of the spinor of particle 2 (index as); it acts on \E' from the 
right; this is also the case of products of 72 matrices, which act on from the right in 
the reverse order: 



~ _ ~ 1 
l2^ll2v^ = ^ai/32(7!^7M)/32a2 ' (^aa/S 



2i 



fa = 1,2) 



(2.15) 



The compatibility condition of Eqs. ( p.l3a| )-( |2.13b|) yields the same constraints as in 
Eqs. (|2.4| ) and (|2.5|) (with ^ and V replaced by \l/ and V , respectively). In particular, 
one has: 

V = V{x^,Pl,p^,Ji,72) ■ (2.16) 

Using the operator Hq [Eq. (plTl) l, Eqs. ( ^l3a| )-( ^l3bD can also be rewritten 
single equation, provided the constraint ( |2.8| ) is used: 



:(7i.pi-mi)(-72.p2-"^2)-f^o ^ 



C{p) 



Si{pi)S2{-p2)Ho 



-V^H , (2.17) 



C{P) 



where Si and S2 are the propagators of fermions 1 and 2, respectively. 

The norm for the internal wave function ip takes the form (in the cm. frame): 



J rf^x Tr Ji _ + 47io720^'( 



2 dV 



2Po 



(2.18) 



with satisfying the hermiticity condition 

= 7io720^7io720 



(2.19) 



Equation ( |2.18 ) shows that, for energy independent potentials (in the cm. frame) the 
norm of ip is not positive definite. In order to ensure positivity, it is sufficient that the 
potential V satisfy the inequality 

1, 



-Tr(\/V) < 1 



(2.20) 



In this case one is allowed to make the wave function transformation 



V^V 



-1/2 



(2.21) 



and to reach a representation where the norm for cm. energy independent potentials is 
the free norm. 

In this respect, the parametrization suggested by Crater and Van Alstine for local 
potentials that commute with 71L72L (hence V"^ = V), 



V 



tanh V 



(2.22) 



satisfies condition (|2.20|) and allows one to bring the equations satisfied by \1/ [Eq. ( p.21| )] 
into forms analogous to the Dirac equation, where each particle appears as placed in 
the external potential created by the other particle, the latter potential having the same 
tensor nature as potential V of Eq. ( |2.22| ). 

We shall henceforth adopt parametrization (|2.22|) for potential V and, according to 
Eq. (p.21|) , shall introduce the wave function transformation 



^ = (cosh\/)^' . 



(2.23) 



[For potentials that do not commute with 71L72L, the generalizations of Eqs. (p.22| ) and 
(|2.23|) and of the corresponding norm are presented in Ref . .] The norm of the internal 
part of the new wave function \l/ then becomes (in the cm. frame): 



1 + 4710720-Po 



dV 



gp2 



2Pn 



(2.24) 



[The relationship between \1/ and ip is the same as in Eq. 

In order to obtain the final eigenvalue equation of ip, one has to decompose the latter 
along 2x2 components. We choose, for this decomposition, the basis of the matrices 1, 
7l, 75 and 7^75: 

^ = ^1 + lLtp2 + 75^3 + 7l75^4 ■ (2.25) 

When potential V is a function of (i.e., not an integral operator) it is possible, for 
some general classes of potential, to eliminate through the wave equations ( p.l3a| )-( p.l3b|) 
three of these components in terms of a fourth one. We choose ip^ for the remaining 
component, which is also a surviving component in the nonrelativistic limit. For a general 
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combination of scalar, pseudoscalar and vector potentials, the final eigenvalue equation 
is of the Pauli— Schrodinger type, the radial momentum operator appearing in it through 
the Laplace operator. For this class of potential, V is decomposed along the 7 matrices 
as: 

T T 

V = V, + 715725^3 + l^lMi^V2 + g'^JU, + , (2.26) 



X 



where Vi is the scalar, V3 the pseudoscalar and V2 the timelike vector potential; U4 and 
T4 are the spacelike vector potentials. [The numbering of the indices and the notations 



are related to the particular projectors to which these potentials are attached |29|. 
After making a change of function of the type 



^3 = Ms 



(2.27) 



the expression of A being given in Ref. 
type equation (in the cm. frame): 



one finds the following Pauli— Schrodinger 



mr — TTin 



4 4 4M2 ^ 4P2 ^ 



> 03 = 



(2.28) 



where the dots stand for the various spin dependent or Ti -dependent terms and L is the 
orbital angular momentum operator. [The complete form of the equation can be found 



in Ref. [29].] We recognize that the scalar potential V\ acts as a modification of the 



total mass M of the fermions through the change M Me^^ while {m\ — is kept 
fixed. The timelike vector potential V2 acts as a modification of the cm. total energy 
Pl through the change Pl^^"^ ■, while (p^^ — P2l) = ("^1 " ""^2) ^^P^ fixed. 

The nonmodification of (m^ — vn^ is simply due to the constraint ( p.8| ). The spacelike 
potential f/4 changes the orbital angular momentum operator from L to Le ^ ^^d the 
combination f/4 + T4 of the spacelike potentials changes the radial momentum operator 
from "Pr to 

p,e-2(f/4 + T4) 

(in the classical limit). The pseudoscalar potential appears 
only in spin dependent and /i^-dependent terms. All these effects are what one expects 
from an external field interpretation of the potentials in the Dirac and Klein— Gordon 
equations and therefore provide additional justification for the parametrization ( p.22|) 



and the decomposition of V along the form (|2.26|) 
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Wave equations can also be written down for a system made of one spin-| and one 



spin-0 particles EH]. We shall not, however, consider such a system in the present work. 
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3 Connection with the Bethe— Salpeter equation 
and the scattering amplitude 



In order to find the connection of the wave equations ( p.la]) -( p.lb|) and (|2.13a|) - (|2.13b| ) with 



the Bethe— Salpeter equation, it is natural to project, with an appropriate weight factor, 
four-dimensional quantities such as Green's functions, scattering amplitudes and wave 
functions on the constraint hypersurface ( p.8|) and to iterate the corresponding integral 
equations around that hypersurface. There is an arbitrariness in the choice of the weight 
factor, but it is limited by several additional conditions, such as correct one-body limit, 
absence of spurious singularities, correct 0{a^) effects, interchange symmetry between 
particles 1 and 2, hermiticity of the potential, etc.. We shall first present the general 
method of approach and then shall discuss the question of the choice of the weight factor. 

Let G{P,p,p') be the four-point Green's function (for bosons or fermions) with total 
momentum P and relative momenta p' and p for the ingoing and outgoing particles, 
respectively [Eqs. ( p^.2| )], and from which the total four-momentum conservation factor 
(27r)^5^(P — P') has been amputated. It satisfies the integral equation: 

G = Go + GqKG , (3.1) 

where Go is the free two-particle propagator, 

Go = G1G2 (3-2) 

[Gi and G2 are the propagators of particles 1 and 2, respectively], and K is the Bethe- 
Salpeter kernel |]1], ^, ^. The off- mass shell scattering amplitude T (from which the total 
four-momentum conservation factor has been amputated) is related to K by the relation: 

T = K + KGoT. (3.3) 

In the vicinity of a bound state pole, G behaves as 

r n(TD '\ ■ (t>{P,P)4>{P,P') p2 J +1 for bosons, 

lim G{P,p,p) = 11] , s = P , r] = < (3.4) 

s-ss s-SB + te 1^ _1 for fermions, 

where (j) is the internal part of the Bethe— Salpeter wave function and is defined from 
hermitian conjugation and antichronological product. 
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We next define the left-projected Green's function on the hypersurface ( |2.8| ): 

G{p,p') = -2t7r6{C{p))G^^G{p.p') , (3.5) 

where Gq^ is the weight factor to be determined below. [When not necessary, total 
four-momentum will not be explicitly written in the arguments.] Notice that G is still 
four-dimensional with respect to its right argument. 

It is now possible to iterate G around G and obtain the integral equation for G. We 
shall use the following definitions and abbreviated notations: 



Sc = SiC{p)) 



g = Go + 2m5cGQcGQ , 



g = +2m5cG^^Go 



fi'O — GqqGq 



C{P) 



(3.6) 



Using in the right-hand side of Eq. (p.5|) the integral equation for G, one obtains for 
the difference G — G: 

G-G = g{l + KG) , (3.7) 

which yields: 

G = {l-gK)-\G + g) . (3.8) 

One then successively uses Eqs. (|3?T|) and (|3.8|) in the right-hand side of Eq. (|3.5| ) to 
obtain: 



G 



-9 



(3.9) 



[l-Kgr' + Kil-gK)-'G 

The second term in the right-hand side of this equation can also be expressed in terms of 
the scattering amplitude with the use of Eq. 



G 



-9 



'l-Kgy' + T{l-gT)-'G 



(3.10) 



We define the constraint theory wave function ijj with the same projection as in Eq. 



^= -2tn5cG,^^, (3.11) 
which leads to the following behavior of G in the neighberhood of a bound state pole: 

ipip^ip') 



lim G(p,p) 

s^sb 



17]- 



S — Sb + I'G 

[t] is defined in Eq. ( ^.4| )], which, in turn, yields the wave equation satisfied by ip: 

^ = -gT{l-gTr'i, . 



(3.12) 



(3.13) 
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We shall also define the reduced wave function ip and amplitude T as: 



T{p,p') 



2P, 



-Tip,p'] 



C{p),C(p') 



The wave equation satisfid by if) takes the form: 

g^^ = -f{l~gof)-'i; 



(3.14a) 
(3.14b) 

(3.15) 



where now the integrations between T, go and are three-dimensional and concern the 
transverse variables [Eqs. ( |2.3[ )], after the constraint ( |2.^ ) has been used for the longitu- 
dinal momenta. 

It is also possible to reconstitute the Green's function G and the Bethe— Salpeter 
wave function (p from the knowledge of the constraint theory Green's function G and 
wave function ip. To this end, we consider Eq. ( p.7| ), in the right-hand side of which we 



replace G by its expression (|3.8|), and use Eq. ( p.9|) for G; we find: 

G = Go + GoT(l - gTy\G + g) , (3.16) 
<P = GoT{l-gT)-'tP . (3.17) 

Notice that in the left argument of the first T in the above equations the constraint ( p.8|) 
is not used. We also dispose for G of Eq. ( |3.8| ), which yields for the relation 

= {l-gXy'tP , (3.18) 

but these equations are more singular than Eqs. (|3.16| ) and ( |3.17D above, and should only 
be used inside integrals. 

It is natural at this stage to define the potential V (also denoted by V for the bosonic 
case) by the Lippmann— Schwinger type relation: 

V = f{l-gof)-' . (3.19) 

Equation ( |3.15|) then becomes: 

g^'i; = -Vi^. (3.20) 

This equation has the same structure as the constraint theory wave equations ( |2.1a|) - 
(|2.1b ) and ( p. 17 ). Comparison of ^q'^ with the corresponding kinematic operators then 
fixes the expression of Gq^ [Eqs. (|3.5| ) and (|3^ )]: 

Ho . (3.21) 



C(P) 
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Let us now examine the hermiticity property of potential V, which should partly be 
guaranteed by the on-mass shell elastic unitarity property of the scattering amplitude, 
through relation (|3.19|) . After a little algebra |jTO|, one finds for the bosonic case the 



following condition on Qq: 

go - To = 2t7rS{bl{s) + p'^)e{s - (mi + m^f) , (3.22) 
where we have defined 

^li.^) = ^(-5 - ("^1 + ^2Y){S - [nil - 7712^) 

= £-l(m; + ,„l) + M^. (3.23) 



The argument of the 5-function above is nothing but the operator Hq defined in Eq. ( p. 11 



which means that solution ( p.21|) for Gqq automatically satisfies the hermiticity-unitarity 
condition (|3.22| ) in the physical region of the amplitude. [Goc should be defined with the 



prescription {Hq + ie) ^.] 

For the fermionic case, the hermiticity-unitarity condition reads: 

% - liLl2L9hiLl2L = 2i7rS{bl{s) +p^'^)0{s - (mi + mg)^) 

x(7i-Pi + "^i)(-72-P2 + "^2) , (3.24) 

which is also satisfied by solution ( p.21| ). 

Other solutions to Eqs. ( p.22 ) or ( p.24 ), different from that given in Eq. ( p.21 



can 



be constructed as well. In principle, different choices of ^0 should not affect the physical 
predictions of the theory, but rather would amount to organizing in different ways the 
resolution of the bound state problem and the summation of the perturbation series, 
putting emphasis on specific approximations. In the peresent work, we stick to the choice 



(|3.21|) , which has the two main advantages of simplicity and of the fact that it reduces 
in x-space to a second-order differential operator. It is only in this case that can we 
device local approximation schemes for potential V and remain in a framework analogous 
to that of nonrelativistic quantum mechanics. Choice (|3.21|) has also the advantages of 
hermiticity, interchange symmetry between particles 1 and 2 and correct one-body limit. 

We shall, however, slightly generalize solution (|3.21|) by allowing in it the presence of 
a finite multiplicative renormalization constant. Such a generalization is necessary in the 
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off-mass shell formalism we are using to avoid the presence of spurious 0(a^)-terms. It 
amounts to requiring that the 1 / r-terms of the potential result solely from the contribution 
of the one-photon exchange diagram. For this purpose we define 

g'o = l~9o , (3.25) 



where go corresponds to solution ( |3.21D and 7 is a constant, expressible as a series of the 



coupling constant. Equation ( |3.2CI| ) then becomes, after replacement of go by g^: 

= -^Vi>, (3.26) 

with 

V = f(l-7^oT)-^ . (3.27) 

The presence of the factor 7 in Eqs. ( p.26| ) and ( |3.27|) leaves enough fredom for accom- 
plishing the desired cancellation. For the simplicity of notation, we shall no longer write 
explicitly the factor 7, but when necessary shall mention its presence. 

The integral equation satisfied by the Green's function G also allows one to construct 



with standard techniques P, O] the scalar product of states. 
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4 Perturbation theory calculation of the potential 



Our main task in the remaining part of this work is the calculation, in perturbation 
theory and within a definite approximation scheme, of potential V from formula ( p.l9| ). 
The perturbation series of V contains, in addition to the usual Feynman diagrams, other 
types of diagram, arising from the presence of the constraint factor qq; we shall call 
these diagrams "constraint diagrams" ; they are obtained from the usual box-ladder type 
diagrams by the replacement of (at least) one pair of fermion and antifermion (or boson) 
propagators by the constraint propagator g [Eq. ( p.6|) ]. The role of these diagrams will be 
to cancel the spurious singularities (referred to the potential) arising from the calculation 
of the amplitude and to yield a potential that has most of the required physical properties. 

It was already shown in the case of fermions |^ that, when one of the particles 
becomes infinitely massive, the contributions of the constraint diagrams cancel those of 
the ladder and crossed ladder diagrams, leaving only the contribution of the one-photon 
exchange diagram. Thus the wave equation of the remaining particle reduces in this limit 
to the Dirac equation in the presence of the static potential created by the infinitely 
massive particle. 

In the present work we limit ourselves to the evaluation of the ladder and crossed ladder 
diagrams in their leading order, neglecting radiative corrections. In this approximation 
the scattering amphtude T [Eq. ( p.l4b| )] can be decomposed as: 



T = ^t(") , (4.1) 

n=l 

where T*^"-* is the partial amplitude corresponding to n exchanged photons. 
Iteration of Eq. (|3.19|) yields for the potential the expansion: 

oo oo 

V ^ = Tj2{goTy, (4.2) 

n=l p=0 

where V^"'^ is that part of the potential which comes from n exchanged photons. The first 
two terms of V are: 

^(1) = fii) ^ (4.3) 

y(2) ^ f{2) ^ fiDg^fW (4.4) 

The corresponding diagrams are represented, for the case of fermions, in Figs. 1 and 2, 
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respectively. [The ingoing particles have now momenta pi and p2 and the outgoing ones 
p[ and The constraint diagram is represented by a box with a cross in it.] 

Our evaluation of the importance of various terms hinges on the infra-red counting 
rules of the bound state problem in QED. Thus, if a is the dimensionless coupling constant 
squared, r [Eq. (|2^ )] counts as 0{a~^), the transverse momentum as 0(a), the 
transferred momentum squared t as 0(q;^), etc.. Furthermore, we shall evaluate terms 
with respect to x-space; therefore, expressions computed in momentum space will undergo 
three-dimensional Fourier transformations and will receive additional O(a^) contributions; 
thus a/t will be counted as 0(q;^), a^/ (— t) as 0{a^), etc.. 

The one-photon exchange contribution, V^^\ is of order and, in three-dimensional 
x-space, is proportional to a/r. As we shall show in Sees. 5 and 6, the two-photon 
exchange contributions, corresponding to V^'^^ [Eq. ( |4.4|) ], yield as a dominant effect an 
0{a'^) term, proportional, in x-space to o? jr^ . It is therefore expected, generalizing these 
results, that the n-photon exchange contributions yield as a dominant effect an 0(q;^") 
term, proportional to (a/r)". We shall show, in Sees. 7,8 and 9, within an approximation 
scheme, that this is indeed the case. Therefore, the sum of leading terms of each V^'-"^ in 
Eq. (f4.2|) will result in a local function of a/r. 

With the above counting rules, there are still nonleading terms, in V^*^") («>!), say, 
that might compete with leading terms of with n! > n. Such terms have, however, a 
different structure than a simple power of a/r; they are essentially momentum dependent 
and cannot be grouped within a function of a/r; therefore, they will not be considered in 
the summation process, although, in a quantitative perturbation theory calculation their 
effects, which begin at O(a^), should be taken into account, like those of the radiative 
corrections. 

To have a consistency check of our summation procedure, we make, in the following 
sections, parallel calculations of the scalar potential (with scalar massless photons) in the 
bosonic and fermionic cases, where the Feynman diagrams have rather different structures. 
One however expects that the spin independent classical part of the potential in the 
fermionic case, when brought in an appropriate representation, should coincide with the 
bosonic potential. This coincidence does occur and provides additional justification for 
our summation rules. 
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5 Two scalar photon exchange contribution in 
bosonic system 



We consider in this section a system of two spin-0 particles, with masses mi and m2, 
respectively, interacting with scalar (massless) photons. The corresponding lagrangian 
density is: 

C = \d^Ad^A + d^<^\d^,^^ - <^\{m\ + 2m^gA + g'A^)<^^ + (1^2). (5.1) 

We have kept analogy with QED and introduced the interactions with the photon with 
the substitutions m\ m\ + gA and m2 m2 + gA. We define: 

Furthermore, we ignore other types of interaction like 0^ or A'^ which arise from renormal- 
ization. We assume that the corresponding renormalized coupling constants have been 
put equal to zero; more generally, these coupling constants are of order g"^ ~ and the 
corresponding Feynman diagrams yield contact potentials that are of order and hence 
can be neglected in comparison to O(a^) effects. 
Our definitions of the propagators are: 

^1(^1) = , G2{P2) = , D{k) = • (5.3) 

In the process (1) + (2) — > (1') + (2') we designate the corresponding momenta by pi, 
P2, p'l and P2, respectively, and introduce the usual momentum variables: 

S = iPl+P2f = {Pl+P2f = = PL 

q = Pi - p'l = P2- P2 , t = , u = {pi- pa)^ . (5.4) 
When constraint (p.8|) is imposed on the external particles, one obtains: 

PiL = p'lL > P2L = P2L , qi = , q^ = q^^ , (5.5) 



with piL and p2L given by Eqs. 

The Feynman diagrams will be calculated with the external particles considered off 
the mass shell, with longitudinal momenta fixed by the bound state mass Pl, through 



relations (2^), while the transverse momenta squared p and p have, according to our 
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counting rules, orders of magnitude of o? . The off- mass shell deviations are represented 
by the quantities 

\2 22 22 \/2 2 /2 2 /2 Ir. r\ 

which are positive and are of order a^. The mean values of p^"^ and p'^^ in the bound 
state are equal and hence we have a similar equality for and A'^: 

<p'r'^> = < ^ ^ < a2 > = < A'2 > . (5.7) 

To the accuracy of the present calculations, where we are evaluating only leading effects 
of a set of Feynman diagrams with a fixed number of photons, we shall not distinguish 
between p^"^ and p''^'^, and between A^ and A'^, and shall use the equalities ( |5.7| ) without 
the mean values. 

We also introduce the off- mass shell analogue of the quantity bl{s) [Eq. (p.23|) ]: 



b\s) = I - lipl+pl) + ^^-jf^ = -p^' , (5.8) 

the second equality arising from the use of Eqs. ( |2.4| ), (p.8|), and (|2.9| ) and of the definitions 
Pi = PiL + P"^^) pI = P2L + P^^- also easily shows the relations: 

^0 = pIl -ml = pIl -ml = -X' + . (5.9) 

Notice that 6^ > 0, while 69 < 0. Keeping only the leading order terms in a, one also has 
the mean value equality: 

< > = 2<b^ > , (5.10) 
which in turn implies the leading order equality: 

(p'^ + p'y = -2(A^ + ^)- (5-11) 

We first calculate the one-photon exchange contribution (Fig. 1). Because of the 
constraint conditions [Eqs. (|5.5|) ], the momentum carried by the photon is transverse, 
Q = {Ql = 0, g^) [q^'^ = — in the cm. frame]; the passage to a;-space is done with 
the three-dimensional Fourier transformation with respect to q^. Taking into account the 
definition of T [Eq. ( p.l4b| )], one finds for the potential V^^^ [Eq. ([4.3|)] (the potential in 
the bosonic case is not tilded): 

t 
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or, m x-space: 

= _^!!^« . (5.13) 
Pl r 

We next calculate the two-photon exchange contributions. In addition to the diagrams 
of Fig. 2, one also has the two diagrams of Fig. 3. We first introduce the notations for 
the various integrals we shall meet during the calculations. They are similar to those 
introduced by Brown and Feynman |^ and Redhead [^] (but our volume element d'^k 



does not contain a {2tt) ^ factor.) We use the abbreviated notations: 

(0) = e + ze , (3) = {q-kf + ie , 

(1) = (pi - kf -ml + ie , (2) = (p2 + kf -mj + ie , 
{-!') = {-p[ - kf -ml + ie , (-2') = (-p'2 + ^? -ml + ie 



(5.14) 



(0)(3)(1)(2) 



(3)(1)(2) ' J (0)(1)(2) 

= / d'k^^^h^ , = / d'k- 



mm) ' J mm) 



(5.15) 



For the crossed diagrams, where (2) is replaced by (—2'), we use the abbreviated 
notations: 



For the integrals resulting from constraint diagrams and according to the definitions 
(|3.21|) , ( p.l4b|) and ( p.l9|) , we have to suppress one of the boson propagator denom- 



inator, (1) or (2), and replace it by 2i7r(5(C)/(2Pi), the operation being symmetric with 
respect to the interchange (1) ^ (2): 



c 2m r il.kl.klkl) 
'^''^^ " 2Pj ^^^^^^ (0)(3)(1) ' 

^^.^i|/^^M(..)MJ. (5.17) 



21 



Details of the calculations of the above integrals are presented in Appendices A and 
B (see also Ref. [^). For the integrals of the diagrams of Fig. 2 we find: 



J 



S 
+ 



^[yX^ 



(l-7)t] 



St 



PI 



ln(- 



(A 



2\2 



s 



PI 



St 
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J(l,-2') 



+0(a'^ In a" 



(5.18) 
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jc = __ 
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(5.19) 
(5.20) 



^3 ri dj /A^ 

10 7[7A2 ~ (1 ~ l)t] \is 

By taking the sum of J, J(l, —2') and J*", one finds that J*^ cancels the dominant 
singularity in J, while the next-to-leading singularity of J cancels the leading singularity 
of J(l, —2'). One obtains: 



J+ J(l,-2') + J' 



c 



2m' 



ln( 



v\l- 



+ ©(a-^ In a" 



(5.21) 



Taking into account the coupling constants and the other multiplicative factors, the 
potential V"*^^-* resulting from the above contribution becomes: 



2 2 4 

m|m2 g 



771 

-^ln(^) + 0(«^lna-^) . 

1712) 1712 



(5.22) 



Pit 7r2(mf 

This is a spurious 0(Q;^)-term, since the potential is known, from other semirelativistic 
calculations, not to possess such a term. It is a two-body effect and cannot be removed by 
one-particle renormalizations. The origin of this term lies in the off-mass shell extrapola- 
tion of the Lippmann— Schwinger formula ( p.l9| ) with the use of Eq. ( p.21| ). We actually 
also calculated the contribution of the three diagrams of Fig. 2 on the mass shell, giving to 
the photon a small mass. We found that the sum of the three terms of Eq. ( ^.21| ) vanishes 
(up to 0(a;'^ In a;~^)) and hence the term ( |5.22|) is absent from the potential. If we admit 
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that at leading orders on-mass shell and off-mass shell calculations should coincide as far 
as the potential is concerned, we find another justification for removing the above term by 
an appropriate renormalization of the Lippmann— Schwinger formula. This prescription 
amounts to requiring that the 1 / r-terms in the potential arise solely from the one-photon 
exchange diagram. We described the corresponding procedure through Eqs. (|3.25| ) to 
( |3.27| ). Choosing to this order 

g"^ mim2 

7 = l+7i, 7i = TT^r^ 2^H — ), 5.23 

27r^ (mf — 7712) 7712 



one can cancel, with the aid of V^^^ [Eq. j fJ^ ], the term \/(2)J [Eq_ (|]22D]: 



V^(2)J + ^^v^(i) = . (5.24) 

Therefore, the net effect of the three diagrams of Fig. 2 can be considered as zero to 
order in the potential V^'^\ 

We next evaluate the contributions of the diagrams of Fig. 3. The corresponding 
integrals are: 

4 

C(^) = + 0(aMna-i) (a = 1,2), (5.25) 

2paLV-t 



and their contribution to the potential V^"^^ is [we recall that there is a symmetry factor 
2 at the four-line vertex]: 



8Pl PlL P2L V-t 



In x-space, the potential V^'^^ takes the form: 



= ^(^ + , (5.27) 

Pl PlL P2L 



and, together with V^(^) [Eq. ( p.l3| )], the potential becomes to order 



V = y« + F(^) = + 1 (!^ + !d)^ . (5.28) 

Pl r Pl Pil P2L 

In summary, the constraint diagram has cancelled the dominant spurious logarithmic 
singularity of the direct diagram, while the remaining subdominant logarithmic singular- 
ities have been cancelled among the direct and crossed diagrams. Among the two types 
of diagram we met in Figs. 2 and 3, it is only the diagrams of the type of Fig. 3 (the 
"seagull" diagrams) that have contributed to the 0{a^) effects. 
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Similar calculations can also be repeated with scalar QED (vector photons). We 
did such calculations in an arbitrary covariant linear gauge (characterized by a parame- 
ter The vector nature of the interaction renders the calculations more complicated. 
Furthermore, the constraint diagram has an ultra-violet divergence and necessitates an 
appropriate subtraction. The final qualitative results remain, however, the same as above: 
the spurious singularities disappear and the leading term is of order a^, with the correct 



coefficient as expected from scalar QED spectroscopy [1T0| , |T6[ and the Breit hamiltonian 
for bosons. We shall not present the details of these calculations in this article and shall 
be content with the relatively simple case of scalar interactions for the bosonic case. 
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6 Scalar and vector interactions with fermions 



We turn in this section to the calculation of the potentials resulting from one- and two- 
sacalar or vector photon exchanges between a fermion with mass mi and an antifermion 
with mass The interaction lagrangian densities are 



SI 



gij^A^i + (1-^2) 



C 



VI 



e^i7/,A^^i + (1 



(6.1) 



respectively, with the common notation for the coupling constant squared, Eq. (|5.2|) and 



a 



The fermion propagators are: 



Si{pi) 



Air 



S2{-P2) 



(6.2) 



(6.3) 



respectively; the scalar photon propagator is given in Eqs. ( |5.3| ), while the vector photon 
propagator is, in the linear gauge characterized by a parameter ^: 



fc2 _|_ 

The constraint factor go [Eqs. (pj), (|3), and (|;T1|)] is: 



go = Si{pi)S2{-p2)Hc 



Cip) 



(6.4) 



(6.5) 



The potential V^^'' that results from one scalar photon exchange is: 

,2 



(1) 



9 



2Prt 



(6.6) 



[t = g-^^, when constraint C is used on the external particles, Eqs. (|5.5| )], which becomes 
in x-space: 



(1) 



a 



s 



2PLr 



For the vector interaction one obtains for V^^^: 

^2 



(1) 



V 



2Pr 



T T 

7i72 (^M^ -^^^) 



t ' t 



(6.7) 



(6.8) 



which becomes in x-space: 



C T T c 

7l 72 (^/.^ - ^uz. o + 



2 2^ ■ 



(6.9) 
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Upon comparing this expression with Eq. (^]^) [notice that the difference between V 
and V [Eq. ( p.22| )] shows up only at third formal order in a] one identifies the vector 
potentials V2, t/4 and T4 in lowest order: 

In the Feynman gauge = 0) one observes the simple relations: 

V2 = Ui, T4 = . (6.11) 

The difference between the timelike vector potential V2 and the spacelike vector po- 
tentials 1/4^ and T4 is that, at each given formal order in a, the latter, because of the 
multiplicative 7^-matrices, contribute to terms that are 0(q;^) times less important than 
the contribution of the former. This feature makes the computation of the potentials U4 
and T4 rather tricky, since they contribute only to nonleading terms, while V2 contributes 
to the leading ones. In the present work we shall only calculate V2. However, if one makes 
the assumption that relations (|6.11|) obtained in lowest order in the Feynman gauge re- 



main valid in higher orders, the knowledge of V2 will allow the reconstitution of the whole 
vector potential in this gauge. 

For the calculation of V^"^^ we need to calculate the contributions of the three diagrams 
of Fig. 2. To this end, we transform the fermion propagators by bringing the 7-matrices 
in the numerators and giving the denominators a bosonic form. Notice that the passage 
from the direct to the crossed diagram is obtained by the replacement p2 + k p'2 — k, 
which, in the bosonic denominators is equivalent to {—P2 + k). We then obtain for the 
scalar potential the following decomposition in terms of the integrals ( |5.15|) - (|5.17|) : 



(7^7^] = ili-Pi + ^i)(-72.P2 + m2)(J + J(l, -2') + J^) 



2Pl (27r)^ 



-(71 .pi + mi)72^(J^ - ^/.(l, -2') + J^') 
-(-72.P2 + m2H{J^ + J,il, -2') + J^) 



+l^Y2{J,u-J,u{l,-2') + Jp 

(7i.pi+mi)J(l,-2')-7i-^(l,-2') 



72.g . (6.12) 



The transverse matrices 7^ are counted as being of order a, since they connect nondom- 
inant components of the wave function to the dominant one. We summarize here the 
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leading contributions in x-space of the integrals of Eq. ( |6.12| ) (see Appendices A and B): 
J + J{l,-2') + = 0{aHna-^) , 

Jl = ^(GW - , Ml, -2') = + , 

72^^.(7^- J^(l, -2') + J^) = 0{aHna-'), 
7r.(J^ + J^(l,-2') + J^) = 0{aHna-'), 

(7i.pi + mi)J(l,-2')-7i--f(l,-2') 72.g = ©(aMna"^) . 

(6.13) 

The first relation above is understood in the sense of Eq. ( ^.24| ). The expressions of G'(i) 
and G^^-* are given in Eqs. (5.25). 



After replacing the matrices ■jn and 72L by their eigenvalues +1 and —1, respec- 
tively, with respect to the dominant component of the wave function and making the 
approximations mi ~ pn, ~ one finally obtains: 



^ 16P? ^ ' ^ ' 



which, in x-space, is equivalent to: 



The scalar potential is then, up to two-photon exchanges: 



2 



For the vector potential one has to add the 7-matrices at the vertices. For the direct 
and constraint diagrams one obtains expressions having the following structure: 71.72 
[ ] 7i-72, where the quantity [ ] is the corresponding expression obtained in the scalar 
interaction case. For the crossed diagram the structure is: 7ia72 [ ] 7i/372 > with [ ] 
representing the corresponding expression of the scalar case. The interchange of the 
matrices 72 and 72 brings in new terms, which, however, turn out to be of order o? In 
compared to the right-hand side of Eq. ( |6.12| ). Therefore, to order in the potential, 
the contribution of the three diagrams can be obtained from the right-hand side of Eq. 
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( |6.12| ) with the global substitution 71.72 [ ] 71.72 and the cancellations met in the scalar 
case operate again. At leading order, one finds for the timelike vector potential: 

or, in x-space: 

2 

Vy^ = -7iL72L^ . (6.18) 

The timelike vector potential is then, up to two-photon exchanges: 

,,2 



Vy = Vi"^ + V^^ = m72.(^ - ^) . (6.19) 



It can be shown that when the two-photon exchange contribution ( |6.18 ) is added 
to the complete one-photon exchange contribution ( |6.1UD -( p.llD and the nonrelativistic 
limit of the wave equations ( p.l3a] )-( p.l3b| ) is taken to order 1/c^, and after a wave function 
transformation that is equivalent to the passage from the Feynman gauge to the Coulomb 



gauge, one obtains the Breit hamiltonian |Q which is known to yield the correct bound 
state spectrum to order a^. This is also an indication that no other 0(a;^)-terms should 
arise from higher order diagrams. (Verification of correct 0(a'^)-terms with Todorov's 



electromagnetic potential was presented in Ref. [45 



Result (|6.15|) was also shown [^] to contain the correct retardation effects, calculated 
in Refs. ^ by differnt methods. 

In summary, the same type of cancellations as for the bosonic case has occurred for 
the fermions, despite the fact that the tensor structure of the Feynman diagrams is more 
complicated. The terms that contribute to the leading [0(q;^)] effects are those coming 
from the longitudinal components of the vector integrals J^. In these integrals, the linear 
contribution of the longitudinal component of the photon momentum of the numerator 
cancels the dominant contribution of one of the bosonic denominators of the fermion 
propagators and transforms the whole integral into a G-type integral ( 5.15| ). All happens 
as if the fermions had effective interactions of the type of Fig. 3. 
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7 Cancellations with constraint diagrams 

The fact that the two-photon exchange diagrams globally produce in leading order (a^), 
both for the bosonic and fermionic cases, local potentials with the correct coefficients, as 
compared to the bound state spectrum analysis to order 1/c^, suggests that such a result 
might also hold with n-photon exchange (n > 2) diagrams; in this case, we should expect 
to have a local potential proportional to (a/r)". It is, however, not possible to show such a 
result in a rigorous way, because of the complexity of the n-photon exchange diagrams. For 
this reason, it is necessary to device an approximation scheme that takes into account the 
properties of the bound state we are considering, as well as the cancellation mechanisms 
present in the two-photon exchange cases. 

The kinematic domain of the bound state is characterized by the fact that the trans- 
verse momenta of the external particles are damped, by an 0(a)-factor, with respect to 
the longitudinal momenta and the masses. Also, the factor of the two-photon ex- 
change potentials represents in momentum space the three-dimensional convolution (up 
to multiplicative coefficients) of the two photon propagators, in which the longitudinal 
momentum squared, kj^, has been neglected; this means that the entire g^-dependence 
of the potential (at leading order) comes from the photons and not from the fermions. 
Finally, the expressions of the potentials we obtained are independent of the off-mass 
shell condition imposed on the bosons and the fermions (i.e., do not depend on pf — ml 
or P2 — ^11,2). 

It is then natural to assume that, in quantities that are expected to be globally free of 
infra-red singularities, the above properties also reflect themselves in the internal boson 
or fermion propagators. We shall therefore consider the approximation that consists in 
linearizing with respect to the internal photon momenta the inverse boson propagators 
and in neglecting in them the momentum transfer, as well as off-mass shell terms. One 
then obtains for the boson propagators the expressions: 

Giipi -k)- ^ \ , ■ , G2{p2 + k) 



2pi.k + it 2p2.k + ie 

G,{-p[ -k)^ I , G2{-p'2 + k) \ . (7.1) 

2pi.k + ze —2p2.k + ie 



This approximation can be considered as a variant of the eikonal approximation |33 



26, p5| adapted to the bound state problem. Here, one expects for the components of k 
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the orders of magnitude 

\kL\ = 0{a) . (7.2) 

[Actually, will come out of the order of 0{a^).] If the propagators ( [7.1|) were multiplied 
by regular terms, we could immediately neglect in them .k'^ , but in many cases the 
multiplicative terms are singular and such an approximation would be misleading before 
the isolation of the singularities. Finally, because of the on-mass shell treatment of the 



boson or fermion propagators in approximations ( |7. 1| ) , the photon should be given a small 
mass to prevent infra-red divergences from occurring at intermediary stages. 

It can be checked by direct calculation that approximations (|7. 1|) yield the correct 



results ( ^.27[ ), ( |6.15D and ( |6.18D . This can also be checked from the results that we shall 



obtain for the general case of n-photon exchanges. 

We shall show in this section that all types of diagram that have constraint diagram 
counterparts are globally cancelled by the latter in leading order of our approximations. 
The main formula we use is a generalization of an identity already used in the eikonal 
approximation (cf. Ref. |2^, Appendix, and Ref. |2^, pp. 116-117). Let (ci, C2, ■ ■ ■ , c„+i) 
be a set of {n+1) numbers; it can be devided into {n+1) independent subsets of n numbers, 
where in each subset one of the Cj's (i = 1, ■ ■ ■ , n -|- 1) is missing. We have the identity: 

n+1 /n+1 \ , /n+1 \ 

E E [ici)-\c[ + 4)-^ . . . (c; + 4 + . . . + c'j~']^ = E ^0 / n , (r.s) 

i=l perm ' \i=l / ' \j=l / 

where (4, 4, ■■■,4) is a permutation of the subset of n numbers (ci, C2,- ■ ■,Cj_i,Cj+i, 

• ■ -yCn+l)- 

Let us consider in the bosonic case an {n + l)-photon exchange diagram of the type of 
Fig. 4, which generalizes those of Fig. 2. We denote by fci, ■ ■ ■ , kn+i the momenta carried 
by the photons; for definiteness we number them according to their departure point on the 
line of boson 1, from the left to the right. To take into account momentum conservation, 
the factor {2TT)^6'^{J2i=i ~ l) must be included in the corresponding integral. Since 

= 0, we have YJi=i = 0; but according to the approximations ( |7.1| ) we may also 
neglect in the boson propagators; hence, we have there the approximation 

n+1 

E = • (7-4) 

i=l 

The total number of the diagrams considered above is (ri -|- 1)!. In a given diagram we 
have n propagators of boson 1 and n propagators of boson 2; only n photon momenta 
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appear on each of the boson hnes (but not the same in general). The (n + 1)! diagrams 
can be divided into (n + 1) sets, where in each set we have n\ permutations of boson 2 



propagators containing the same set of n photon momenta. With approximations ( |7.1| ), 
the sum of all these propagators has the structure of the left-hand side of Eq. (|7.3|) , with 
Cj = {2p2-ki + ie); hence, we obtain the factor 



'2,n+l 



n+1 

Y^{2p2.ki + ie) 

i=l 

n+1 

Xl{2p2.kj + ie) 



(7.5) 



\_Y{2p2-kj + It) 



Using Eq. ( [7.4| ), we can rewrite /2,n+i in the form: 

n 

h,n+i = ^ I ; — + 

E 

(7.6) 

Repetition of the above operation on /2,„, /2,n-i, etc., yields the final result: 

/2,n+i = {-2^nr5{2p2.k,)S{2p2.{k^ + k^)) ■ ■ ■ 5{2p2. k,) . (7.7) 

i=l 

The 5-functions can then be used, upon integrations on the /cjl's, in the boson 1 
propagators. Denoting by the total contribution of the above diagrams (without 

the coupling constants and other multiplicative coefficients), we obtain: 

1 1 1 

X 



[-2pT.k1 + ze) {-2pT.{kT + kl) + le) {-2pT . kj + ze) ' 

(7.8) 

where /x is a small mass given to the photon, to avoid infra-red divergence. Also, the 
5-functions ( [7.7| ) yield for the kfj^^-teims orders of magnitude of O(a^), which accounts for 
their omission in front of 0(a^)-terms in approximations ( [7.1|) and in photon propagators. 

On the other hand, we should also take into account the contributions of constraint 
diagrams that are associated with the diagrams considered above. From the general 
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formula ( [4. 21 ) we obtain for the potential corresponding to n + 1 photon exchanges the 
expansion: 

n 

yin+l) ^ f(n+l) f^'-^^gof^'-'^go ■ . . f^^p)gQf^^p+^'^ , n>l , n>l , 

p=l riH hrp+i=n+l 

(7.9) 

where, in the generic term of the sum, the constraint factor appears p times. A typical 
diagram, where p = 2, is shown in Fig. 5. The constraint diagrams corresponding to the 
diagrams of the type of Fig. 4 are those for which the amplitudes T'^^^^ in Eq. ( [7.9| ) are also 
of the same type, with exchanged photons. Now the longitudinal momentum transfer 
between the external particles and a constraint factor or between two successive constraint 
factors is zero; furthermore, according to approximations ( |7. 1| ) , we can also neglect trans- 
verse momentum transfers in the bosonic propagators between external particles and 
constraint factors and between two successive constraint factors. Hence, for the partial 
amplitudes T*^^*) (r^ > 2) in Eq. (|7.9| ) we obtain integrals J'-'"*-' that have the same structure 
as in Eq. ( [7. 8] ), with n replaced by rj — 1, except for the factor (27r)^5^ which is relative 
to the whole diagram. Next, the constraint factor go provides the coefficient (2111) / {2Pi) 
multiplied by the boson propagator contribution l/(— 2p-^. X^jL^^^ ^'^^ + ie), and so 
forth. Therefore, the constraint diagrams yield globally the same expressions as 
[Eq. (|7.8|) ], except for a sign factor equal to (— 1)*^ if % occurs p times in the expansion 
(|7.9| ). Moreover, the p factors ^0 niay appear in {ip^ different configurations. The total 
contribution of diagrams of the type of Fig. 4 and of their associated constraint diagrams 
is then proportional to 



j(n+l) 



i + E(-ir(?) 



(1 - 1)"+^ = , n > 1 . (7.10) 

p=i 

This result generalizes the one obtained in the two-photon exchange case. 

With diagrams involving more than two-photon exchanges, we can meet configurations 
that are mixtures of diagrams of the types of Fig. 4 and of generalizations of Fig. 3. In 
this case the permutational procedure on a given line should be applied to those parts 
of the diagrams for which the momentum configuration on the other line is unchanged. 
Furthermore, to show the cancellation mechanism by constraint diagrams, not all propa- 
gators need be permuted; it is sufficient to focus attention on those propagators that do 
not correspond to the basis of a seagull diagram (of the type of Fig. 3), since it is only 
these that may have constraint diagram counterparts. A simple example illustrates this 
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situation. Let us consider the six diagrams of Fig. 6. We apply on line 2 the permuta- 
tional method for the first two diagrams [(a) and (b)]. The result is a (5-function which 
yields the constraint diagram (a) of Fig. 7 with an opposite sign and hence a mutual 
cancellation occurs. Similarly, graphs (c) and (d) of Fig. 6 are cancelled by the constraint 
graph (b) of Fig. 7. 

The application of the permutations on line 2 in diagrams (e) and (f) of Fig. 6 
yields also a 5-function, but the resulting expression does not have a constraint diagram 
counterpart. It is proportional to: 

Sip2-k2) r ■■ TT^T^ • (7-11) 

After an integration with respect to k2L this expression becomes proportional to: 
1/1 1 \ 



(7.12) 



p'^.k^ \-piLkiL - p^-kl + ie -piLhi - p'^-kj - PLP'^.k^/p2L + ie^ 
Finally, by making in the second term the changes of variable 

kiL ^ k[L = k,L + Pl^-^ , k^L ^ k',^ = k,L - Pl^-^ , (7.13) 

PilP2L PilP2L 

the whole integral becomes proportional to: 
1 1 



p'^./cj {-piikiL - p'^.kf + i 



D{kl, + kDDikl, + kD - D{k[l + kDDik'i, + kf 



(7.14) 

With the use of assumptions ( |7.2| ), the difference of the two terms inside the brackets 
yields a non-leading term that can be neglected. Therefore, the six diagrams of Fig. 6 
are cancelled, at leading order, by the two constraint diagrams of Fig. 7. 

This result can be generalized to more complicated diagrams. Use of the permutational 
procedure on a given line produces a certain number of ^-functions, but only a smaller 
number of them survives at leading order; these are precisely those which appear in the 
constraint diagrams and are cancelled by them. 

We thus end up with the conclusion that the only diagrams that survive (at leading 
order) for the calculation of the potential in the bosonic case are those which do not have 
constraint diagram counterparts. These are generalizations of the diagrams of Fig. 3, 
some examples of which are presented in Figs. 8 and 9. 
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The above study can also be applied to the fermionic case. For the fermion propagators 
we use the approximations: 

Si{pi - h) ~ — — ; — — [(71LP1L + mi) - 7iLfciL] , 
—2pi.ki + le 

S2{-{P2 + h)) ^ 7^ ; — -— [{-12LP2L + m2) - ^21^1] ■ (7.15) 

Zp2.K2 + 

In the numerators we have neglected the transverse momenta and kf (but not in the 
denominators); in products of two fermion propagators we neglect in the numerator the 
quadratic terms k1 (of the same ki)- It can be checked that approximations ( [7.15| ) yield 
the correct results ( |6.15| ) and (|6.18|) of the two-photon exchange case. 

With fermions, we have only diagrams of the type of Fig. 4, but the presence of the 
kiLS in the numerators of their propagators gives rise to a situation that is very similar 
to that of the bosons. The products of (n — l) pairs of fermion propagators (for n-photon 
diagrams) can be decomposed into three types of terms. The first type of term does not 
contain any kn in the numerator. It behaves as in the bosonic case with diagrams of 
the type of Fig. 4 and leads to an expression that is proportional to J*^"^ [Eq. ([7. 81)] and 
is cancelled by the corresponding constraint diagram contributions. In the second type 
of terms some of the kiis (but not all) appear in the numerator. The analysis of these 
terms is similar to that of the mixtures of diagrams of the types of Fig. 3 and Fig. 4 we 
considered in the bosonic case (but with the additional possibility of having more than 
two photons at the effective seagull vertices) and their contributions are again cancelled by 
corresponding contributions of constraint diagrams. In the third type of terms products 
of (n — l) independent combinations of the kn's appear in the numerator. These do 
not have constraint diagram counterparts, since there kn = (i = 1, ■ ■ ■ , n). Therefore, 
these terms are the only surviving parts of the above diagrams for the calculation of the 
potential. 
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8 The potential in the bosonic case 



We calculate in this section the scalar potential for two-boson systems. According to the 
analysis of Sec. 7 the diagrams that contribute in higher orders are those of the types 
of Figs. 8 and 9. The diagrams of Fig. 8 correspond to an odd number of exchanged 
photons, while those of Fig. 9 to an even number of exchanged photons. 

We begin by considering diagrams of the type of Fig. 8. Let (2n + 1), with n > 
1, be the number of exchanged photons. The total number of such diagrams is {n + 
l)!(n + 1)!. We denote by ki,k2,---, k2n+i the momenta carried by the photons, and 
for definiteness we number them starting from the photon that leaves the boson line 
1, and then follow continuously the other photons. A global momentum conservation 
factor, {2TT)'^6^(J2^2:t^ ki — g), must also be put in the integrals corresponding to these 
diagrams. The different diagrams involve permutations of the boson momenta, which 
concern, however, independent sets on line 1 and line 2, respectively. Thus, according 
to our numbering of photon lines, the set of internal momenta appearing on line 1 is 
{ki,k2 + ^3,- ■ ■,k2n + k2n+i), whilc the set of internal momenta appearing on line 2 is 
{ki + ^2,^3 + k^,- ■ ■,k2n-i + k2n,k2n+i)- Furthermore, each set satisfies the equation: 

2n+l 

Y^k^L = 0. (8.1) 

i=l 

Since the integrals we are now calculating are of the regular type, we can, from the 
start, neglect the transverse momenta kf in the boson propagators (however, this is not 
compulsory for the following calculations). Application of the permutation formulas ( [7.3|) 
and ( |7.5|) - (|7.7|) for the longitudinal momenta independently on lines 1 and 2 yields for 
the boson propagators the result (notice that formulas ( [7.5| )-( [r7|) applied on line 1 with 
propagators of the type (— 2j9i.A;- + ie) still produce (— 2i7r)-factors): 

1 1 

/£5 = (^2in?"-——-——T]6(k,L), (8.2) 
2n+i I ; (2pii)" (2p2l)" l\ 

and the corresponding potential becomes: 

Y{2n+1) ^ ^(_2imi^)(-2zm2^)(-2i^2)2n.2n.2n+l^_.)2n^l 



/2n+l j3 UT 1 \ 2n+l 

(n^^)(2.mi:*j-.^). (8.3) 
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= -^7^ i-X"" . (8.4) 



yielding in x-space: 

..,'>r,^^^ 4___ 

2Pl {PilP2lY \rJ 

Notice that this formula is also applicable to the case of one-photon exchange [Eq. (|5.13|) ], 
for which it provides the exact result. 

We next consider the diagrams of Fig. 9, with {2n + 2) exchanged photons {n > 0). 
Here, we have to distinguish between the two types of diagram, according to whether the 
photon lines begin and end on line 1 (Fig. 9a) or on line 2 (Fig.9b). We designate by 
y{i)(2n+2) g^^^ ^(2){2n+2) corresponding potentials, respectively. Let us, for definiteness, 
consider the diagrams of the first type (Fig. 9a). The number of such diagrams is 
|(n + 2)!(n + 1)!, where the factor 1/2 comes from the existing symmetry with respect to 
fixed positions of starting and ending photon lines. The remaining part of the calculation 
parallels that of {2n + 1) exchanged photons. We have on line 1 the set of {n + 2) 
independent internal momenta {ki,k2 + k^,- ■ ■,fc2„ + ^2^+1, ^2^+2) and on line 2 the set 
of {n + 1) independent momenta {ki + ^2,^3 + k^,- ■ ■,/c2n+i + ^2^+2)5 each satisfying Eq. 
(|8T| ) (with 2n + 1 replaced by 2r2 + 2). Since the permutational procedure necessitates 
(n + 2)! permutations on line 1 and (n + 1)! permutations on line 2, we have to put the 
combinatorial factor 1/2 in front of the final result. We find: 

Vimn+2) ^ ^(_2,^^^)2(_22(72)2n+1^2„+1^2n+2^_^^2n+l 1 1 1 



2Pl' ^ ^ ' ^ ' 2 (2pii)"+i (2p2l) 

\ (27r)3(fcP+Z6), 



X 



^2n+2 j3 UT i \ 2n+2 



1.5) 



yielding in x-space: 



y(l){2n+2) 



4Pl V \/PiLP2Lr J 



p^{2)(2n+2) jg obtained from Eq. ( P^ with the exchanges 1^2. 

Collecting all the contributions calculated above, we obtain for the total potential: 



^ _ 2mim2 ^ ^ 

\/PlLP2Li / 

2 |2\/ \ 2 nr, / \ 2n 



Pl r „=o \VP^LP2Lr) 



Pl WPiLP2Lr) ^0 \ ^PiLP2Lr) 
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which on summation becomes: 

^ 2'mim2 a ^ f 'm1p2L + 'm'2PiL 2mim2a\ f a / ( ^ / 



Pl r \ Pl Pl r J \ ^/pEpoIr J / \ ^/PiW2Lr J 

(8.8) 

This expression can be further simphfied. We used for the evaluation of higher order 
diagrams approximations where we neglected (in the regular expressions) the transverse 
momenta of the external particles (as they contribute to nonleading terms). This is 
equivalent to using the static limit in these diagrams. Therefore we are entitled to use the 
equations of motion in the static limit in the corresponding expressions to eliminate energy 
factors in terms of the masses and energy independent potentials. Using the equations of 
motion (|2.1a| )- (|2.1b|) in the static limit [p^"^ = 0) with V given by Eq. (|8.8|) , we can, after 



a few algebraic manipulations, reexpress the second piece of V (which involves higher 
order terms only) in terms of simpler quantities. The final result for V is: 



2mim2 a a 



2 



V = (8.9) 

Pt. r 



Potential V above has the same structure as the scalar potential proposed by Crater 



and Van Alstine ||2^ on the basis of an extension of the minimal substitution rules intro- 
duced by Todorov who identified the two-particle motion in the cm. frame to that 
of a fictitious particle with reduced mass and energy defined, respectively, as: 



where 6g [Eq. ( ^.23] )] is the mass-shell invariant relative momentum squared, correspond- 



ing to the cm. total energy ^/s = Pl. It can be checked, after rewriting the Klein— Gordon 
operator Hq [Eq. ( p.ll| )] in the form 



Ho = E^-ml, + p^^ , (8.11) 
that potential (|8.9|) is generated in the wave equations ( p.la|) - (p.lb| ) by the substitution: 

ml. ^ {mw + Sf , (8.12) 

with S = —a/r. 
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9 The potential in the fermionic case 



We turn in this section to the calculation of the scalar and timelike vector potentials 
relative to fermion-antifermion systems. The fermion propagators are considered in ap- 
proximations ( [7.15|) . According to the analysis, at the end of Sec. 7, of n-photon exchange 
diagrams, the only surviving terms, at leading order, are those containing in the numer- 
ator products of (n — 1) independent combinations of the longitudinal photon momenta 
kiL- These arise from pairs of fermion and antifermion propagators containing in their 
numerator one ki^. The permutational procedure already used in other instances pro- 
duces (5-functions which allow each kn to cancel the denominator of the corresponding 
fermion propagator that did not participate in the permutational operation. Hence, there 
remains at the end the convolution of n photon propagators. The combinatorial analysis 
is much simplified if we use these cancellations (with the appropriate coefficients) prior 
to the permutational procedure. 

In a diagram with n exchanged photons (Fig. 4) we have several possibilities to choose 
the (n — l) independent combinations of ki^s. Considering for definiteness line 1, we can 
keep the {n — 1) k^ factors in the numerators of the propagators; then we have to remove 
the (n — l) kiL factors from the numerators of the propagators of line 2. We can also 
keep (n — 2) kn factors on line 1 and keep one kn on line 2, the choice of the latter being 
completely fixed by the complementarity condition to the {n — 2) ki^s of line 1: the set of 
the {n — 1) fcji's globally kept must always form an independent set of vectors. In general, 
we have the possibility of keeping p kiLS on line 1 and the complementary [n — 1 — p) 
kiis on line 2. 

Let then p be the number of ki^s kept on line 1. This choice can be made in p ''") 
different ways; however, no freedom is left for the choice of the set of (n — l—p) complemen- 
tary kiis on line 2. During the permutational procedure, we need {p+l)\ permutations to 
obtain p (5-functions on line 2 from the propagators not having ki^s in their numerator; 
similarly, we need {n — p)\ permutations to produce {n — p — 1) 5-functions on line 1. 
Taking into account the total number of diagrams with n exchanged photons, which is n!, 
we obtain the combinatorial factor relative to the above permutational procedure: 

where F{a,b; c; z) is the hypergeometric function |4^, the value of which, for z = 1, is 
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known in terms of F-functions. 

After making the approximations {'JilPil+^i) — '^PiL and (— 72LP2L+'Ti2) — 2p2L, the 
integrations on the ki^s produce the global factor ((— 2i7r)/(2PL))"^^. The corresponding 
potential becomes, in the scalar interaction case: 

(n) _ ^ / ■^.\2n-2(n-l)-n ""^ \ " ^ (2^^)' 

1 



^(_,^)2n,2(n-l)^ 

2Pi^ V2Pl 



X 



n 



in + l)!n! 



yielding in x-space: 



(n) 



L (27r)3 (fcP + 
(2n)! 



(9.2) 



(9.3) 



^ (n + l)!n! V 2Pir. 

This expression is also valid for the one-photon exchange case [Eq. ( |6.?| )], for which it 



provides the exact result. The two-photon exchange result [Eq. ( |6.15|) is also reproduced. 
The total scalar potential becomes: 



(2n)! / a 



n=l 



{n + l)\n\ \2PLr 



2a 



1 + 



2a 
Plt 



(9.4) 



It satisfies the inequality (|2.20|) , and hence parametrizations ( p.22|) and (|2.26| ) can be 
used. One obtains: 



Vs = tanh Vi 
1 



- In 1 + -— 

4 V Plt 



(9.5) 



The timelike vector potential is obtained from the scalar one, by multiplying, for an n- 
photon exchange diagram, t/J"^ by (— 1)"(7il72l)", the factor —1 coming from the vector 
photon propagator and the factor 7il72L from the vertices. Replacing {n — 1) factors 
7il72L by their eigenvalue —1 at leading order, we end up with the global factor — 7il72l 
with respect to V^'^K We thus find: 

(2n)! 



a 



{n + l)\n\ V 2Pir 
leading for the total potential to the expression: 

'l 



llLl2L 



(9.6) 



V 



-llLl2L 



1 + ^ 
^ ^ Plt 



1 + 



1 + ^ 



(9.7) 
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In terms of parametrizations ( |2.22| ) and ( ^.261 ) we have: 



Vy = tanh(7ii72LV2) , 

V, = ]ln(l + ^] . (9i 



4 V Plt. 

If we adopt for the electromagnetic potential the hypothesis that in the Feynman 
gauge relations ( |6.11D , obtained in lowest order, remain also valid in higher orders, we can 
reconstitute the entire potential: 

V, = U, = ^ln(l + |^) , = 0, 

Vv = tanh(7i.72V2) • (9.9) 
This potential is the fermionic generalization of Todorov's potential [|lOl, obtained in 



two spin-0 particle systems from Eqs. (|8.10| ) and (|8.11|) with the substitution rule 



E' ^ {E + VoY , (9.10) 

where Vq = a/r. Indeed, if we go back to Eq. (|2.28|) , which essentially represents the 
squared two-body Dirac equation, and replace there the potentials V2, f/4 and T4 by their 
expressions ( |9.9D , we find that in the classical part of this equation the electromagnetic 
interaction is introduced with the substitution rule ( |9.1(JD . 



As to the scalar potential, it does not correspond in its form (p.4|)-(p75D to the substi- 



tution rule ( ^.121 ), but can be made compatible with it if in the higher order terms the 
classical static equations of motion are used, as in the bosonic case [Eq. ( |8.9|) ]. To this 
end, let us replace Vi by its expression (|9.5| ) in Eq. ( p.28| ) and retain in the latter the 
classical static (and spin independent) part {p^ = 0): 

4 + 4P2 4 e ~ ' ^ ^ ^ 

from which we deduce the relation 

Pl ~ Me^^i , (9.12) 



to be used in higher order terms. Equation ( p.ll[ ) can also be rewritten in the form 



^e^y^{l-e-^V^Y - 0. (9.13) 
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Using in the last term, which concerns higher order terms than the first, Eq. ( |9.2| ) for the 
factor M^e^^^i we find: 

where 5" is defined as: 

S = . (9.15) 

r 



The interaction dependent part in Eq. ( p.l4|) has the same form as in the bosonic case 



[Eq. (|8.9|) ], after a similar procedure was used there, and corresponds to the substitution 
rule ( WT% ). 

This agreement can be considered as a consistency check for our approximation scheme, 
applied to bosons and fermions with different diagrammatic structures. As one would 
naturally expect, the classical parts of the corresponding effective interactions do coincide. 

Adopting Eq. (|9.14|) as the final form of the classical part of the effective interaction. 



one can reconstruct the new expression of Vi through the identification of the classical 
static part of Eq. (^^) with Eq. (|9A|). One finds: 



2 , 2mim2 c , c2\^ , f 2 , 2mim2 2 



mi + + S'j + [mi + + S 



(9.16) 



with S defined in Eq. (|9.15| ). 

The results obtained thus far can be generalized to include the case of a mixture 
of scalar and vector interactions. The calculations can be repeated as above, with the 
difference that every scalar photon propagator (of the scalar interaction case) is now 
replaced by an effective propagator that is the sum of the scalar and vector propagators 
(including the couplings at the vertices): 

g^Ds ^ g^Ds + e^iiLl2LDv • (9.17) 

It is this sum which is factorized in the Feynman diagram integrals, while the fermion 
propagator part remains unaffected. Concerning the determination of the tensor nature 
of the interference terms, the approximations utilized in the fermionic case (in particular, 
the replacements in several places of the matrix products 7il72l by —1) do not leave 
enough predictivity about this question. It turns out, however, that the analysis of the 
same problem in the bosonic case is more predictive and suggestive of the representation 
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of the interference terms by effective scalar interactions. We shall adopt this property 
also for the fermionic case. 

When the summation is done with the effective propagator (|9.17| ), the pure vector 
potential, which has the form (|9.7|) - (|9.8|) , can be isolated from the total potential and 
7il72L replaced in the rest by —1. The result is: 



V = tanh{Vi + 7il72l^2 
V2 



a 
r 



S 



Pl 



a 
r 



a 



3_ 

4:71 



(9.18) 



Including then the spacelike potentials U4 and T4 with the hypothesis ( p.9|) , we obtain: 



V = tanh(T/i + 7iL72L"t4 + 7i -72 U4) 
U, = V2 = ^ln(^l + -^K) 



(9.19) 



We again can transform the scalar potential with the aid of the classical static equa- 
tions of motion, as in Eqs. ( 9.11| )-( 9.16| ). We end up with the following expression of 
Fi: 



1 



In 



1 



MJI + -irVo 



Pl 



P, 



Pl 



^1 2n , 2 , 2mim2 2 
+ { "^2(1 + -^^0) + S + S 



a 



Vn = - 



a 



S = 



(9.20) 



with V2 and U4 given by Eqs. (|9.19|) . The classical parts of the above potentials satisfy 
Todorov's minimal substitution rules fig, defined in Eqs. (glOD-(pID, (gJ^) and ( CT )- 
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10 Reconstitution of the Bethe— Salpeter wave 
function 



In this section we shall reconstitute, in the fermionic case, the Bethe-Salpeter wave func- 
tion in the framework of the local approximation that we used throughout this work. The 
exact reconstitution formula is given by Eq. ( p.l7| ), where in the left argument of the first 
T the constraint ( |2.8| ) is not used. A first approximation consists in replacing the factor 
T(l — gT)^^ by potential V [Eq. (|3.19|) ] in which, however, pi is no longer submitted to 



the constraint (|2.8|) . This is easily done by leaving, in the expression V{t,PL), calculated 
in the local approximation, the part q1 of t{= q\ + g-^^) as a free variable, not fixed at its 



zero value. More precisely, qi = Vl — v'l— Vl 



7712) /{2Pl), where p'^ acts on the 



right on ip and hence is submitted to the constraint ( p.8|) . We therefore obtain: 

2Pr, rd^p"" ' '-^ -2 



^i.P) = Sl{pi)S2{-p2)- 



p 



(10.1) 



A further approximation can be imposed by neglecting altogether the g| dependence 
of V in Eq. (|10.1|) . In this case the wave equation ( ^.201) , together with Eq. (|6.5| ), can be 
used and one obtains: 

2Pl 



-S1S2 



1 -1 



SiS2H(i 



C{P) 



^ ■ 



:io.2) 



The relative longitudinal momentum is now contained only in the two propagators 5*1 and 
Formulas ( |10.1| ) or ( |10.2| ) can be used in the calculations of transition amplitudes, 
form factors and decay coupling constants. 

Of particular interest is the expression of the Bethe— Salpeter wave function at the 
zero value of the relative longitudinal coordinate [it enters in the calculation of decay 
coupling constants]: 



{xL = 0) 



dpL 

271 



4>{pl 



(10.3) 



Adopting approximation ( 10. 2| ), integration on pi concerns only the propagators Si and 
5*2. One finds: 



dpL 
2tt 



Si{pi)S{-p2] 



APlHo l^^^-^^ + -i)(-72-P2 + m2) + 



El E2' 



72l(7i-Pi + "^1) + 7iL (-72.^2 + m2) 



43 



where Ei and E2 are defined as 

E2 = 



Cip) 



(10.4) 



El = \/mi- 



;io.5) 



and constraint ( p.8|) is applied in the right-hand side [pn and p2L are given by Eqs. (|2.9|) ]. 

In the framework of local potential approximation, we can replace the nonlocal op- 
erators El and E2 by the eigenvalues pn and p2L, respectively [Eqs. ( p. 91) ]. Equations 
(|TU3)-(|TU3) then yield: 



{XL = 0,X') = \1 



2VpiL P2lJ 



72l(-72-P2 - rn2) + 7il(7i.Pi - mi) 



+7iL72L7r(7i-Pi - "^i)(-72-P2 - ^2) \ ^/'(x'^) . 
-no ) 

After the use of the wave equations ( p.l3a| )-( ^.13b|) , Eq. ( p.0.6| ) becomes: 



(10.6) 



0(xl = 0, x^) = 1 + 7il72L^ 





( 1 


-^) 






P2l) 



72l(7i-Pi +mi) 



+7il(-72-P2 + ^2) 



(10.7) 



Finally, using parametrization ( p.22| ) and the wave function transformation ( p.23|) one 
obtains: 



(^[xl = 0, x'^) 



2\piL P2l) 



72l(7i-Pi + rni] 



+7il(-72-P2 + ^2) 



sinh y| ip{x^) , 



(10. 



where ip is normalized according to Eq. ( p.24|) and V is of the type (|2.26|) and is a function 
of x^ and P^. 
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11 Conclusion 



Constraint theory, applied to two-particle systems, provides a natural basis for a mani- 
festly covariant three-dimensional reduction of the Bethe— Salpeter equation. The two- 
body potential is related, as in the quasipotential approach, to the off-mass shell scattering 
amplitude by means of a Lippmann— Schwinger type equation and is calculable in terms 
of Feynman diagrams. Perturbation theory is reorganized with the presence of "constraint 
diagrams" that appear in the course of the three-dimensional reduction process. 

The two-photon exchange diagrams (for scalar interactions or vector interactions con- 
sidered in the Feynman gauge) are globally free of spurious infrared singularities and yield 
at leading order local potentials proportional to (a/r)^, where r is the three-dimensional 
relative distance in the cm. frame; these, together with the one-photon exchange con- 
tribution, produce the correct bound state spectra to order a'^. The n-photon exchange 
diagram contributions were evaluated in an approximation scheme that is a variant of the 
eikonal approximation adapted to the bound state problem; the latter produce at leading 
order local potentials proportional to (a/r)". The series of leading order potentials were 
summed, producing total potentials that are functions of r and of the cm. total energy. 
The potentials calculated for vector and scalar interactions are equivalent, for their classi- 
cal parts, to those obtained by Todorov and Crater and Van Alstine |^ on the basis 
of mimimal substitution rules. 

In summing the series of leading order potentials, ambiguities might arise from possible 
incorporation or omission of nonleading terms. A crucial restriction imposed on such 
ambiguities comes from an inequality [( p.20| )] that should satisfy the fermionic potential 
in order to ensure the positivity of the norm of the wave function. The potentials obtained 
in this work do satisfy the above inequality and hence do not lead to spurious violations 
of the positivity of the norm. 

Another question that arises concerns the utility of the potential formalism itself for the 
evaluation of the bound state spectrum: would not the scattering amplitude, calculated in 
the eikonal approximation, provide the bound state spectrum through the determination 
of the positions of its poles? The accuracy of the evaluation of the potential in the present 
work goes beyond that of the conventional eikonal approximation. For the bosonic case for 
instance, the latter sums diagrams of the type of Fig. 4. From the potential theory view- 
point, the net effect of the contribution of such diagrams is equivalent to the one-photon 
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exchange contribution (Fig. 1). It is the summation of the chains of seagull diagrams of 
the types of Figs. 3,8 and 9, not taken into account in conventional eikonal approximation, 
that provides the genuine multiphoton-exchange contributions to the potential. Similarly, 
in the fermionic case, the contributions isolated in the present work for the evaluation 
of the potential are not considered in conventional eikonal approximation; futhermore, 
the latter would produce a spin independent bound state spectrum. On the other hand, 
the potential evaluated above is reinjected into relativistic wave equations that take into 
account more accurately spin dependent effects. We conclude that potential formalism is 
necessary for a detailed probe of the relativistic bound state spectrum. 

The generality of the summation method, in which the photon propagators, except for 
the counting rules, do not play any active role, leaves open the possibility of extending the 
above calculations to other types of interaction, like those corresponding to the exchanges 
of massive particles or those corresponding to effective confining interactions. Another 
domain of application might concern the study of the gauge dependence problem of the 
electromagnetic two-body potential. 

The fact that the constraint theory wave equations can be reduced to a single Pauli- 
Schrodinger type equation |2^, provides the complementary basis for their simple appli- 
cability to relativistic bound state problems. 
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A Calculation of the scalar integrals 



In this appendix we calculate the scalar integrals defined in Eqs. ( ^.ISD - dS.lTI ). They have 
already been evaluated in the literature |^2|, but generally on the mass shell, with a 
small mass given to the photon. Since we are using an off-mass shell formalism, we present 
here some details of the calculations. Definitions of variables and approximations were 
introduced at the beginning of Sec. 5, while the counting rules of the orders of magnitude 
were explained in Sec. 4; we recall that these are referred to x-space. We use for the 
calculation of the integrals the Feynman parametrization method. 

We first consider the integral J. After an integration on one of the parameters, we 
find: 



J 



1 rl 



JO 



d'-yda 



[7A2-(l-7)t] 



A2 + 7[ap? + (1 - a)pl - a{l - a)s] 



(A.l) 



We next integrate with respect to cr; the result is an arctan function, from which we 
isolate the dominat part by expressing it with respect to the inverse of its argument. In 
the nondominant part we make the change of variable x = 



IT 

7S 



^. We find: 



J 



J + J 



(A.2) 



with 



?7r 



3 .1 



and 



J 



dx 



In 



1 - Jl 



{A2)2 \ 



bit 



1 + A 1 



St J^s [x^-{{X-f-hlt)/{st)] V-^W£) + arctan(J: 



(A2)2 

(A.3) 



(A.4) 



where 



- ^. 



(A.5) 



The integration domain of x in Eq. ( |A.4| ) is separated into two parts by a point lying 
in the region of the x's of the order of a^/^. Appropriate approximations of the integrand 
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in each of the domains lead to the following expression of J: 



J 



St 



1 1 



+ 0{a^ Ina"^; 



(A.6) 



J(l, —2') is obtained from Eq. ( |A.1| ) by the replacement of s by -u [we are working in 
the approximations P2 = p\ and A'^ = A^]. One first integrates with respect to 7: 

r2 ri da 



J(l,-2') = - — 



t Jo 



apl + (1 - a)pl - a(l - a)u + A2(A2 + t)/t 



(A.7) 



In the bound state domain we are considering, we have the following approximate expres- 
sions for s and u: 



Pi + h^ + ^pI + 



I Pi + ^ Pi] 1 + 



PiVl' 



u 



pIpI 



(A. 



which allow us to write: 



apl + (1 - a)pl - a(l - a)M ~ ( ^pl + (T(^pf - v^pl 

^1 + ^2? 



+a{l-a) 



PiPl 



(A.9) 



the second term in the right-hand side is negligible in front of the first and this leads to 
the following expression of J(l, —2'): 
2z7r2 



J(l,-2') 



St 



2V/5i /^a/ V St 



+ 0(aMna-i: 



(A.IO) 



The constraint integral J*", which is three-dimensional, is also calculated with the 
Feynman parametrization method. After one integration one finds: 



J' 



c 



it: 



3 rl 



-J 



(A.ll) 



48 



The structure of the F-integrals is very similar to that of the J-integrals ( |A.1| ) and 
the same types of approximations are apphed. We obtain: 



/ 



f + F, 

s . 

s 



2m^ 







7s s 



S ^ 

— arctan \ / — ^ , 
¥ V % 



(A.12) 
(A.13) 



n-^ ( 2 ^ uii \ arctan — + arctan — , 



+ 0(aMna-^) , 
(A.14) 



-2') = -ztt' 



[o-pi + (1 - (7)p2 - - cr)^i] 



In 



' (yp\ + (1 — cr)p2 — cr(l — (y)u 



A2 



r_i/^ Ml 

2VA + /52y ""^s^ (A-/52)V/3i 



^ ^^lnA--^ln/52 

P2 



+0(a^lna-^) 



= — , 

s Jo /A^ ^ 



ITT' 



3 H 



(i7 



(A.15) 
(A.16) 



Notice that with the approximations pf = pf, p\ = p'2, = A'^ we have used, the 
F-integrals are independent of t and hence 



H 



H(l.-2') = Ffl,-2' 



5 y 5 



H 



c 



(A.17) 



The G-integrals can be expressed in terms of Euler's dilogarithm or the Spence func- 
tion. One finds (without the approximations pf = pf and pi = P2): 



it: 



^-2t{pl+pf)+t^ + {pl-pfY 
the dominant part of it being: 



+ 0(a^lna-^) 



(A.18) 



2piL\/-t 

G^^^ is obtained from the above equations by the replacements 1^2. 

G^^^ does not have a crossed-diagram counterpart, since it depends only on the par- 
ticle 1 momenta. On the other hand, the constraint integral G^ [Eq. ( p. 17] )] does not 
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correspond to an explicit diagram in the expansion ( [4.4| ), but rather may appear through 
integrals involving fermions or vector interactions of bosons. One finds: 

= /• , , - ^ . (A.20) 



We emphasize that the instantaneous approximation in the photon propagators (set- 
ting there = 0) would not produce the results found above. For instance, for the 
G- integral one would find twice the correct result ( [A.19| ). A detailed analysis of the in- 
tegral shows that the photon propagator poles in ki contribute with a factor —1/2 with 
respect to the contribution of the boson propagator poles. However, the linearization ap- 
proximation of the inverse bosonic propagators, as used starting from Sec. 7, makes the 
instantaneous approximation compatible with the leading order results for the potential. 
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B Calculation of the vector and tensor integrals 



Most of the vector and tensor integrals can be calculated by algebraic decompositions 
with the aid of the vectors pi, p2, q, etc., in terms of the scalar integrals calculated in 
Appendix A, except for those decompositions that might involve infra-red divergences, in 
which case a direct calculation is necessary. We present here the final results [notations 
are those of Appendix A] : 



\-(-- 









FM = 0) = pIh , 
^2 r , 1 

s 



62 



+ 0(aMna-i) 



F^(l, -2') = - — \p[^^^ + -Pim) 



ITT 



+ 
+ 



2/5i/92V/3i-/?2 



Pl(^2 
1 



2/3i/?2V/?i-/32 



1 , , 1 

I +C»(q;^ In q;-^) , 

I + 0(aMna-^) , 



?7r 



q^F +_p^^-(^-arctan. 



q,F^+p';H 
^3 r^/A2 



62 



-|)+0(aMna-i) 



A2 



I2(p2+pf) 2t(p2+p'2) 



P\{p\-l^lf 

2tipl+p?y 



1) 



c 



(B.l) 
(B.2) 



(B.3) 



(B.4) 



(B.5) 



(B.6) 



(B.7) 
(B.8) 
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Glf-* is obtained from G^^^ by the replacements pi —p2, v'\ ^ ^v'l ^"^^ q ^ q: 



For the J^-integrals one finds: 

J, = P,l^,{G^'^ - G^'^) + q^liJ + ]{F - H)] 



(B.9) 



(B.IO) 



J^(l, -2') = -P,^ + + 2(A^ + ^) J(l, -2') - ( -2') + H{1, -T 



{PiL-P2l),.2 , ^ 



(A^ + -) J(l, -2') - _2') 



J(l,-2') + -(F(l,-2')-i/(l,-2') 



(B.ll) 



c 



-{p^ + p' 



(B.12) 



In showing Eqs. ( |6.13|) for the transverse components of the J^-integrals, one uses in 
some places the leading order equality ( p. lip . 

For the J^i,-integrals it is preferable to first make a decomposition along the vector 



integrals already calculated and gJ^J. For the simplicity of presentation, we use certain 



relations previously obtained (such as F = H, etc.) and the leading order equality ( p.ll| ), 
and neglect some nonleading terms. We find: 



^(PilG«+P2lG(2)) + AV 



1 


■ 1 




V2Pl 



1 



'2{pT+p'T\2 



(PilG«+P2lG(2^) + 2^ 



-^(mG«+P2LG(^^)-(2A^ + |)J 



+F + H 



2t 



(-A^ + hj + f-H- ^{pilG^'^ + P2lG^'^) 



(B.13) 
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[/ and H are defined in Eq. ( [B.l| ),] 

1 



p2 



i(G«+G(2)(-2')) + (y + t)J.(l,-2') 



{PlL - P2l) 



"(r + P" -2') + 2^M^^(i> -2') 
+.Jjfi- ^"^"i"^V (i,-20 + Q(c.^), 



2P2 



(B.14) 



2G 



c 



1 



(B.15) 



[H^ is defined in Eq. (|B.5| 
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Figures 



Fig. 1. One-photon exchange diagram. 

Fig. 2. Two-photon exchange diagrams; the "constraint diagram" is denoted by a cross. 
Fig. 3. Two-photon exchange diagrams contributing to the bosonic case. 
Fig. 4. A multiphoton exchange diagram. 

Fig. 5. A typical "constraint diagram" where the constraint factor qq [Eq. ( |4.2| )] appears 
twice (p = 2). 

Fig. 6. Examples of mixtures of elementary diagrams in the bosonic case. 
Fig. 7. "Constraint diagrams" associated with diagrams of Fig. 6. 

Fig. 8. Typical diagram with an odd number of exchanged scalar photons contributing 
to the potential in the bosonic case. 

Fig. 9. Typical diagrams with an even number of exchanged scalar photons contributing 
to the potential in the bosonic case. 



57 




58 




Fig. 5 
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